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Phase Aberrations in Diffraction Microscopy 
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In coherent X-ray diffraction microscopy the diffraction pattern generated by a sample illuminated 
with coherent x-rays is recorded, and a computer algorithm recovers the unmeasured phases to 
synthesize an image. By avoiding the use of a lens the resolution is limited, in principle, only by 
the largest scattering angles recorded. However, the imaging task is shifted from the experiment to 
the computer, and the algorithm's ability to recover meaningful images in the presence of noise and 
limited prior knowledge may produce aberrations in the reconstructed image. We analyze the low 
order aberrations produced by our phase retrieval algorithms. We present two methods to improve 
the accuracy and stability of reconstructions. 
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I. INTRODUCTION 



A new imaging technique has emerged in recent years 
that can overcome many limitations of light, electron, 
and X-ray microscopy. Coherent X-ray Diffraction Mi- 
croscopy (CXDM)i promises to enable the study of thick 
objects at high resolution. In this technique one records 
the 3D diffraction pattern generated by a sample illu- 
minated with coherent x-rays, and as in x-ray crystal- 
lography a computer recovers the unmeasured phases. 
This is done by alternately applying constraints such as 
the measured intensity in reciprocal space and the ob- 
ject support — the region where the object is assumed to 
be different from — in real space. This corresponds to 
defining the envelope of a molecule in crystallography. In 
our implementation the support is periodically updated 
based on the current object estimate^. 

By avoiding the use of a lens, the experimental re- 
quirements are greatly reduced, and the resolution be- 
comes limited only by the radiation damage 3 ^. However 
the imaging task is shifted from the experiment to the 
computer, and the technique may be limited by our un- 
derstanding of the phase recovery process as well as the 
algorithm's ability to recover meaningful images in the 
presence of noise and limited prior knowledge. 

Recently we have presented experimental results of 
high-resolution 3D X-ray diffraction imaging of a well- 
characterized test object to demonstrate the practical 
application of these advances 5 ^. Here we extend the 
analyis of image reconstruction and determine low-order 
phase errors (essentially image aberrations) that can oc- 
cur when reconstructing general complex- valued images. 
We present two methods to improve the accuracy and 
stability of reconstructions. 



II. COHERENT X-RAY DIFFRACTION 




FIG. 1: (a) SEM image of the 3D test object. Scalebar is 1 
micron, (b) Infinite depth-of-focus image reconstructed from 
a central section of the 3D coherent X-ray diffraction data. 

Three-dimensional coherent X-ray diffraction data 
were collected at the Advanced Light Source^ from a 
test object that consisted of 50-nm diameter gold spheres 
located on a 2.5-/im-wide silicon nitride pyramid (Fig. 
GJi). A bare CCD lo,cated in the far field recorded the 
diffraction patterns with a pixel sampling that was more 
than 4 times the Shannon sampling rate for the (phased) 
complex amplitudes. Diffraction patterns were collected 
for many sample orientations over an angular range of 
129 ° . These were interpolated onto a 3D grid. We re- 
constructed a full 3D image by performing phase retrieval 
on the entire 3D diffraction dataset (i.e. the iterations 
involved three-dimensional FFTs). The resulting volume 
image reveals the structure of the object in all three di- 
mensions and can be visualized in many ways including 
projections through the data, slices (tomographs), or iso- 
surface rendering of the data. 

In addition to 3D images, we perform much anal- 
ysis and algorithm development on 2D datasets. For 
the work in this paper we choose central plane sec- 
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tions extracted from the 3D diffraction pattern. By the 
Fourier projection theorem, the image formed from a cen- 
tral section is an infinite depth-of-focus projection im- 
age (Fig. Qd). We carry out ab initio image recon- 
structions using the Relaxed Averaged Alternating Re- 
flections (RAAR) algorithm 9 with the "Shrinkwrap" dy- 
namic support constraint 2 . Details of the algorithm pa- 
rameters used are given in Chapman 5 . 



III. RESOLUTION ANALYSIS 




The phase retrieval process recovers the diffraction 
phases with limited accuracy, due to factors including 
SNR of the diffraction amplitudes, missing data, the in- 
consistency of constraints, and systematic errors in the 
data (such as errors in interpolation). These errors in 
phase reduce the resolution of the synthesized image. 
With a complex image a loose support constraint will 
lead to unconstrained low-order aberrations. As is well 
known an object could be shifted by a few pixels each 
time we reconstruct, which is equivalent to a varying 
linear phase ramp in reciprocal space. In addition to 
this shift low order phase variations, such as defocus and 
astigmatism can also be unconstrained if the aberrated 
object fits inside the support. One way to quantify the 
effect of these phase variations is to determine the vari- 
ation in retrieved phases as a function of resolution 10 
10 -'. Given a reconstructed image g(x) obtained by phase 
retrieval starting from random phases, and its Fourier 
transform G = \G\exp{iip(q)}, we define the phase re- 
trieval transfer function by 



PRTF (q) = |(exp{^(q)})| 



g(q) 

|G(q)| 



(1) 



with (G) the average over the complex diffraction ampli- 
tudes of many reconstructed images starting from ran- 
dom phases. Where the phases are random and com- 
pletely uncorrclatcd, the average will approach zero. 
Thus, the ratio is effectively a transfer function for 
the phase retrieval process, and the average image (the 
Fourier tranform of (G) ) is the best estimate of the im- 
age: spatial frequencies are weighted by the confidence 
in which their phases are known. 

In our case when reconstructing complex 2D images, 
with low frequencies missing due to the beamstop, we 
have observed that phase retrieval from independent ran- 
dom starts may differ by a phase vortex (right or left 
handed), centered at the zero spatial frequency (Fig. 
We find that we can improve the estimate of the image 
by separating out the vortex modest. These phase vor- 
tices are due to stagnation of the phase retrieval process. 
Other phase vortices can appear near local minima of the 
measured intensities, and our method of separating so- 
lutions will fail to detect vortices not centered near the 
beamstop. In order to remove these vortex aberrations 
we modified the reconstruction algorithm as follows: (i) 



FIG. 2: Phase difference between two reconstructions in re- 
ciprocal space, showing a phase vortex between two solutions 
in the far field. The center of the vortex is at q = 0, and the 
half- width of the phase map is q = 0.048 nm _1 . 



Average n independent reconstructions which will likely 
average out the phase vortex modes but will also smooth 
the resulting image, reducing the resolution, (ii) Refine 
this averaged image by inputting it to the RAAR 9 algo- 
rithm and carrying out 200 iterations. Using this "av- 
eraged RAAR" algorithm we reduced the probability of 
recovering an image with phase vortex mode from 40% to 
15%, resulting in an improvement of the PRTF by almost 
a factor of two. 

half period (nm) 
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FIG. 3: Top: Phase retrieval transfer function of the two al- 
gorithms. PRTF=1 represents stable phases. The averaged 
RAAR algorithm significantly improves the stability of the 
reconstructions. Bottom: PRTF after removing the vortices 
centered at q = 0. The averaged RAAR shows marked im- 
provements at high frequencies. 

We compute the final image, and the PRTF, by av- 
eraging 1000 such reconstructions (Fig. 0). Before av- 
eraging many images we make sure that they are not 
shifted with respect to one another by finding the linear 
phase ramp that minimizes the difference between their 
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Fourier transforms. Fluctuations of the linear phase term 
indicate fluctuations in positions. Fluctuations in higher 
order polynomial phase terms indicate that phase aber- 
rations are present in the reconstructions. 

To quantify the instabilities of these low order phase 
modes, we find the low order phase modes (focus, astig- 
matism, coma, up to a polynomial of order n p ) that min- 
imize the difference between each new reconstruction G n 
and the first recovered image Go- This is done by mini- 
mizing 

X = \ Q (q) - Gn (q) exp {ip(q)}\ 2 (2) 
q 

with the 2D polynomial defined by coefficients pij as 



FIG. 4: Demonstration of reference-enhanced diffraction 
imaging. A Pt dot was deposited near the sample (a coccolith) 
with a focused ion beam instrument. The Fourier transform 
of the diffraction intensities (right) can be used to determine 
the support which can constrain the low-order phase aberra- 
tions. 



P (q) = Pi,j 1: 



X 



(3) 



i,j=0 



with q x ^ y = q x .y/2 max (q x .y)- The linear terms repre- 
senting shifts in real space are found using the method 
described by Fienupii, while higher order terms are ob- 
tained by fitting the phase difference, arg(GjG„), to the 
higher order 2D polynomial terms and iterating until the 
correction is less than 1 ° . The fluctuations of the second 
order polynomial coefficients are obtained by calculat- 
ing their standard deviation among 1000 reconstructions, 
and we find that 



Po,o Pl,0 P2A 
*'<l I P0,1 Pl,l P2,l 
P0,2 Pi, 2 P2,2 



2.26 0.31 0.15 
0.53 0.61 
0.14 



(4) 



The linear terms (pi,o Po,i) represent a shift of 
[0.31,0.53]/2tt = [0.049,0.085] pixels in real space corre- 
sponding to 0.5 and 0.8 nm shifts in x and y. The degree 
of defocus phase variation depends on (j>2,o +Po,2)/2, and 
the real-space defocus variation is given by: 



Sz = 



X 



AnNA 2 



std {p 2 ,0 +P0,2} 



(5) 



We have NA=0.084 and A=1.65 nm, giving Sz = 11.3 
nm. Note that this defocus variation represents an in- 
stability of the phase retrieval process and does not cor- 
respond to an optical effect of focusing through a thick 
object. In this case all voxels of the 3D images or pixels 
of the 2D projection images are equally aberrated by this 
effect. 

An additional method to further reduce these instabil- 
ities is to use a small reference point near the specimen. 
During the retrieval process the image of the reference 
point is forced to be small with a tight support. This con- 
strains the aberrations at this image point, and hence at 
all image points. The reference point has the additional 



advantage of providing a hologram of the specimen (Fig. 
|1J which can be used to provide the object support, or 
even the desired image. 

To quantify our ability to recover unmeasured intensi- 
ties (for example behind the beamstop) we use the nor- 
malized standard deviation 

a 2 (q) = (|G(q)-G(q)| 2 )/|(G(q))| 2 . (6) 
We define a transfer function, based on a 2 as: 



TF(q) = 



1 



\A + <r 2 (q) 



l(G(q))| 



G(q) 



(7) 



which has the desired properties that transfer function is 
unity for a 2 = and zero for a 2 — oo. Eqn. (7J reduces 
to the PRTF in the regions of q where |G| is measured. 

An algorithm that always recovers the same phases 
does not necessarily recover the correct ones. Another 
requirement is that the recovered image is constrained 
in the region called support: g (x) — 0, x ^ S. If 
this condition is satisfied the Fourier modulus condition 
(|G| = I 1 / 2 ) is unlikely to be satisfied in the presence 
of noise. We can quantify deviations from the mea- 
sured values by an i?-factor (similar to that used in 
crystallography 12 ) by 



'RF 



(q) 



|G(q)|-VT(q) 



V/(q) 



(8) 



and 
[1 + 



its 



'RF 



(q 



related transfer function RFTF (q) = 

1 12 

, which is plotted in Fig. for a re- 



constructed image. 

IV. CONCLUSION 

We have performed a characterization of high- 
resolution imaging of an isolated 3D object by ab ini- 
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FIG. 5: The R-Factor Transfer Function (RFTF) of a recon- 
structed image, showing excellent agreement with the mea- 
sured diffraction intensities. 

tio phase retrieval of the coherent X-ray diffraction, and 
examined metrics to allow the quality of image recon- 
structions to be assessed. 

The phase retrieval process does not produce unique 
images, in that varying low-order phase modes arise, akin 
to aberrations in an imaging system. Other than the tilt 
terms, the low-order phase aberrations discussed here will 
be reduced in case of a real object (for which only an- 
tisymmetric terms are allowed) and will not be present 
when a real-space positivity constraint can be imposed, 
since defocusing or otherwise aberrating an image causes 
it to be complex. However, in the case of samples consist- 



ing of more than one material (such as biological samples) 
the object cannot be considered positive and we must re- 
duce the effects of aberrations. We have proposed two 
methods of overcoming limitations of computer recon- 
struction: in order to improve the stability of the recon- 
structions we average several reconstructed images and 
use the result to feed a new round of phase retrieval. 
From an experimental point of view, the use of a refer- 
ence point, or other well-defined object, should enable us 
to greatly reduce low order phase aberrations. 
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